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We study scaling behavior of a chiral order parameter performing a simulation of two-flavor QCD 
with improved Wilson quarks. It has been shown that the scaling behavior of the chiral order 
parameter defined by a Ward-Takahashi identity agrees with the scaling function of the three- 
dimensional 0(4) spin model at zero chemical potential. We extend the scaling study to finite 
density QCD. Calculating derivatives of the chiral order parameter with respect to the chemical 
potential in two-flavor QCD, the scaling property of chiral phase transition is discussed in the 
low density region. We moreover calculate the curvature of the phase boundary of the chirl phase 
transition in the temperature and chemical potential plane assuming the 0(4) scaling relation. 
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1. Introduction 

It is important to study the nature of QCD phase transition at high temperature and low density 
for understanding the evolution of the early universe and analyzing experimental data obtained by 
the heavy ion corrosions. In order to extract unambiguous signals for the QCD phase transition 
from the heavy ion collisions, quantitative calculations directly from the first principles of QCD 
are indispensable. At present, the lattice QCD simulation is the only systematic method to do so. 

Many important properties of finite temperature quark matter have been uncovered by lattice 
simulations. Recent developments of computational techniques enabled us to extend the study to 
small chemical potentials. Most lattice QCD studies at finite temperature (T) and chemical poten- 
tial {}x q ) have been performed using staggered-type quark actions with the fourth-root trick of the 
quark determinant. Because the theoretical base for the fourth-root trick is not confirmed, it is in- 
dispensable to carry out simulations adopting different lattice quark actions to control and estimate 
systematic errors due to lattice discretization. Another unsatisfactory point of the staggered-type 
quarks is that, it was very difficult until recently [|l|] to confirm the scaling properties around the 
critical point, which is universal to the three-dimensional 0(4) spin model for 2-flavor QCD, as 
expected from the effective sigma model analysis 1 . This may suggest large lattice artifacts in the 
results with staggered-type quarks. 

Several years ago, the CP-PACS Collaboration has studied finite-temperature QCD using the 
clover-improved Wilson quark action coupled with the RG-improved Iwasaki glue [Q, B]. With 
2-flavors of dynamical quarks, the phase structure, the transition temperature and the equation of 
state have been investigated. In contrast to the case of staggered-type quarks, both with the standard 
Wilson quark action [Q] and with the clover-improved Wilson quark action the subtracted chiral 
condensate shows the scaling behavior with the critical exponents and scaling function of the 0(4) 
spin model in a rather wide range of the parameter space near the chiral phase transition. 

We would like to highlight scaling properties of the QCD phase transition by numerical simu- 
lations with dynamical Wilson quarks. As shown in Fig. |] (left) for \L q = 0, the phase transition is 
expected to be first order when the up and down quark masses (m uc i) and strange quark mass (m s ) 
are sufficiently large or small, and becomes crossover in the intermediate region between them. We 
expect a second order transition in the chiral limit of 2-flavor QCD, i.e. m u d = m q = 0,m s = oo, 
and the nature of the transition changes as m s decreases. It becomes first order below the trictitical 
point (m^). When \x q ^ 0, the phase diagram in the (T,jj, q ) plane is expected as Fig. |l| (right) for 
2-flavor QCD with m q = and m q ^ 0. At low density, the phase transition is expected to be second 
order for m q = and crossover for m q ^ 0. To confirm this phase structure, the scaling study in the 
crossover region is important. Recently, the scaling study with the physical strange mass has been 
done using a staggered-type quark action in Ref. [p, and the scaling property including a chemical 
potential has been also discussed [S]. Therefore, it is worth revisiting the scaling study using the 
Wilson-type quark action. In particular, we want to clarify the scaling property at finite \i q . 

In this report, we study 2-flavor QCD, as a first step. Previous results for 0(4) scaling at 
jl q = are summarized in Sec. |2[ We then extend the discussion to finite density in Sec. ||. The 
scaling properties are tested by numerical simulations in Sec. ||. The conclusion is given in Sec. |5| 

'Precisely speaking, the chiral transition of 2-flavor QCD with staggered-type quarks is universal to the 0(2) spin 
model, whose scaling properties are close to those of the 0(4) model. 
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2 flavor QCD quenched QCD 




Figure 1: Left: Mass-dependence of the nature of (2+l)-fiavor QCD. Right: Expected phase structure of 
2-fiavor QCD at finite temperature and density for m q = and m q ^ 0. 

2. Scaling behavior of the chiral order parameter for 2-flavor QCD 

The chiral phase transition in 2-fiavor QCD at finite temperature is expected to have the same 
critical properties as the 3 -dimensional 0(4) spin model. Hence, the scaling tests of the chiral order 
parameter (yy) provide us with a useful way to study universality properties of the chiral phase 
transition for 2-flavor QCD. 

The order parameter of the 0(4) spin model is given by the magnetization M. In the vicinity 
of the second order critical point, M satisfies the following scaling relation: 

M/h l / s =f(t/h l /y), (2.i) 

where h is the external magnetic field, t is the reduced temperature which is defined by t = 
(T — r c |/ )= o)/7c|/i=o- The critical exponents have the values l/y = l/(j85) = 0.537(7) and 1/5 = 
0.2061(9), and f(x) is the scaling function. In 2-flavor QCD, we identify M ~ {iffy), h ~ 2m q a 
and t ~ j3 — j8 rt . Here, j8 c , is the critical point in the chiral limit, and a is the lattice spacing. We 
compare the scaling functions of 2-flavor QCD and 0(4) spin model. 

A careful treatment must be required because the chiral symmetry is explicitly broken for 
Wilson quarks at finite a. In Ref. ^J, it have been shown that the 0(4) scaling of Eq. ( |2~T| ) 
is well satisfied when one defines the quark mass m q a and the chiral order parameter (yy) by 
Ward-Takahashi identities The quark mass is defined by the following correlation function for 
temporal direction: 

2m q a = -m PS (A»P(0))/ (P(t)P(O)), (2.2) 

where P and are the pseudo-scalar and vector meson operators, respectively, mps is the pseudo- 
scalar meson mass, and the bar means the spatial average. The chiral order parameter is given by 
the following equation: 

2m a a , , . , ... 2m a a(2K) 2 , , , , N , 
(W) = T^rL^W)) = V N 1 (tr (irVsO-Vs)) . (2.3) 
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Figure 2: Left: The ciral order parameter in 2-flavor QCD as a function of the Ward-identity quark mass for 
each j3 . Right: 0(4) scaling plot of the chiral order parameter. 



Here, D is the quark matrix. The quark mass and the chiral condensate satisfy the Ward-Takahashi 
identity in the continuum limit: 

(d^A^(x)P(x)) -2m q a{P{x)P{x')) = <5(x-x')(w)- (2.4) 

We plot the data of the chiral condensate as a function of the quark mass in the left panel 
of Fig. ||. The closed symbols are obtained in the study of Ref. [Q] by CP-PACS Collaboration, 
and the open symbols are new data in this study. (See Sec. |4|) We then reconstruct the data into 
M/h l l s = f{t/h { ly) for J8 = 1.80 - 1.95 and 2m q a = 0.0 - 0.9, which are shown in Fig. | (right). 
The critical exponents of the 0(4) spin model are used. The dashed line is the scaling function 
obtained by the 0(4) spin model in Ref. [J7|]. We adjust three fit parameters in this analysis. One is 
the critical value of j8, which is j8 rt = 1.462(66) in Ref. [Q], and the others are used for adjusting 
the scales of the horizontal and vertical-axes to the scaling function of the 0(4) spin model. This 
scaling plot indicates that the scaling function of 2-flavor QCD is consistent with that of the 0(4) 
spin model. 

3. Scaling property at finite density 

Next, we discuss the universality property of 2-flavor QCD at finite density. It is expected that, 
when \i q is increased from in the chiral limit, the second order phase transition changes to first 
order at the tricritical point, as shown in the right panel of Fig. 1. The 0(4) scaling behavior is 
expected around the second order transition line even at finite densities because the action has the 
chiral symmetry in the chiral limit. At small \i q , we may identify the scaling variables, t and h, by 
the following way for finite density QCD: 

M=(yff ¥ ), t = p-p ct + j(^y, h = 2m q a. (3.1) 
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Table 1: Simulation parameters for m,p$/my = 0.65 (left) and mpg/wv — 0-80 (right) [p 



p 


K 


T/Tpc 


Traj. 


1.80 


0.145127 


1.07(4) 


5000 


1.85 


0.143502 


1.18(4) 


5000 


1.90 


0.141849 


1.32(5) 


5000 


1.95 


0.140472 


1.48(5) 


5000 



j8 


K 


T/Tpc 


Traj. 


1.80 


0.141139 


0.93(5) 


6000 


1.85 


0.140070 


0.99(5) 


6000 


1.90 


0.138817 


1.08(5) 


6000 


1.95 


0.137716 


1.20(6) 


6000 



Here, c is the curvature of the critical line in the (j3,n q /T) plane. On the other hand, h does not 
have a /^-dependent term at low density. Because the critical line is expected to run along the 
m q = axis in the low density region of the {m q ,pL q /T) plane, h is zero at m q = 0. 

It is worth discussing the phase structure investigating the scaling behavior of the chiral order 
parameter. One of the easiest ways for confirming this scaling property is to calculate the second 
derivative of the chiral order parameter. We expect the following scaling properties: 



d 2 M 



d(n q /T) 2 



dM 
dt 



dM/dt 



fcl/S-l/y 



df(x) 



n„=o 



dx 



c=t/h l /y 



(3.2) 



We study the second derivative performing numerical simulations of 2-flavor QCD (Nf = 2) at 
jj, q = 0. The derivative of the chiral order parameter is computed by the following way. We define 



% = (2K) 



2 d n tr(D- l y 5 D- l r5 ) 



d{}X q a) n 



£ n =N f 



d"lndetD 

d{v q a) n ' 



itives of the chiral condensate are given by 





M 9 =0 


2m a a 


d(W) 




2m q a 


d(lh/T) 




nM 


d 2 (yy) 




2m q a 


d(Li q /T) 2 




NfNf 



(3.3) 

(3.4) 

(3.5) 
(3.6) 



2 _ 2^i^ - ^o^ 2 + 2^ 2 ) = (JF 2 - &*sfi) , (3.7) 

where we used the properties that s$ n and & n are zero for odd n's at \i q = 0. The operators, =S„, 
can be calculated by a random noise method. 

4. Numerical simulations for the calculation of the chiral condensate 

We calculate the second derivative of the chiral order parameter for 2-flavor QCD at \i q = 
using the configurations obtained in the simulations of Refs. pi 0]. The RG-improved gauge action 
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Figure 3: Left: 0(4) scaling plot of the second derivative of the chiral order parameter. Right: The curvature 
of the critical line in the (j3 , piq/T) plane determined by the chiral order parameter. 



and the Nf = 2 clover-improved Wilson quark action are adopted. The simulation parameters are 
summarized in Table []} The details of the simulation are given in Ref. [§, ^]. 

We used the random noise method of Ref. \$\ for the calculation of the inverse of the matrices 



in Eq. ( |3.3| ). As we emphasized in Ref. [|9j], it is important to apply the noise method only for the 
space index and to solve the inverse exactly for the spin and color indices without applying the 
noise method to obtain reliable results. We choose the number of noise vectors 50 for each color 
and spin indices. The Ward-identity quark mass obtained by the zero temperature simulation in 
Ref. [Q] is used for the scaling plot. It has been shown in Ref. [Q] that the range j8 > 1.95 is out of 
the scaling region, i.e. the data of {yy) do not agree with the 0(4) scaling function in that range. 
We thus plot the results of the second derivative of the chiral condensate at j8 < 1.95 in the left 
panel of Fig. ||. The colored line is the derivative of the scaling function c{df(x)/dx), which is 
obtained by performing a numerical differentiation of f(x) in the right panel of Fig. || numerically. 
The additional parameter c is chosen as c = 0.02,0.03,0.04 and 0.05 from bottom to top. These 
results are roughly consistent with the expected scaling behavior. Because the error is still large, it 
is too early for drawing a definite conclusion. In particular, the error becomes large as j8 decreases. 
But, the result suggests that the chiral order parameter satisfies the scaling relation in Eq. ( |3.2| ). 

We then calculate the curvature of the critical line assuming the scaling relation. The ratio 
of (d 2 M /d{pL q /T) 2 )h- Y l 8+l ly and df(x)/dx\ x=t/hi/y is plotted in Fig. | (right). This gives c = 
d 2 fi ct /d(pL q /T) 2 , i.e. the curvature of the critical line in the chiral limit. Taking the average, we 
obtain d 2 fi ct /d{lx q /T) 2 = —0.0328(35). Moreover, the curvature of the critical temperature T c (jx q ) 
at \l q = can be calculated by 

1 d 2 T c _ d 2 ^ /dfi 
T c d(ll q /T) 2 d(ll q /T) 2 / da' 

In this calculation, we need a(df}/da) at f} ct in the chiral limit. The value of a(dfi/da) for the 
lattice action which we used have been calculated in Ref. [J3J] for the pseudo-scalar-vector mass 
ratio nips/my > 0.65. If we try to estimate adopting the value near the pseudo-critical temperature 
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at m PS /m v = 0.65, which is a(dp/da) « -0.5, the value of {\/T c )(d 2 T c /d(}l q /T) 2 ) « -0.07. 
This value is smaller than that we expect from the data of the chemical freeze out, but this result 
is similar to the result in Ref. Further studies are, of course, necessary. However, it is found 
that this scaling analysis provides a systematic way to investigate the phase boundary of the finite 
density QCD in the low density region. 

5. Summary 

We studied the scaling property of the chiral phase transition in the low density region of finite 
temperature and density QCD. Assuming the phase structure shown in Fig. [I] (right), the scaling 
function in the low density region is discussed. To confirm the scaling property, we calculated 
the second derivative of the chiral order parameter performing a simulation of 2-flavor QCD with 
improved Wilson quarks and compared the results with an expected scaling function. Because the 
chiral symmetry is explicitly broken for the Wilson quark at finite lattice spacing, we used the chiral 
order parameter defined by a Ward-Takahashi identity for the scaling analysis. The scaling behavior 
of the chiral order parameter for 2-flavor QCD is roughly consistent with the scaling behavior of 
the corresponding spin model even at finite density. 

Moreover, we calculated the curvature of the critical line in the chiral limit assuming the 
scaling relation. It is important to reduce the error by increasing the statistics and the number 
of noise vectors for more precise study. However, we found that this scaling analysis provides a 
systematic way to investigate the phase boundary of finite density QCD in the low density region. 

This work is in part supported by Grants-in-Aid of the Japanese Ministry of Education, Cul- 
ture, Sports, Science and Technology, (Nos. 21340049, 22740168, 22840020, 20340047) and by 
the Grant-in-Aid for Scientific Research on Innovative Areas (Nos. 20105001, 20105003). 

References 

[1] S. Ejiri, F. Karsch, E. Laermann, C. Miao, S. Mukherjee, P. Petreczky, C. Schmidt, W. Soeldner, and 
W. Unger, Phys. Rev. D 80 (2009) 094505 [arXiv =0 90 9.5122]. 

[2] A. Ali Khan et al. (CP-PACS Collaboration), Phys. Rev. D 63 (2000) 034502 [hep- 1 at / 0008011]. 

[3] A. Ali Khan et al. (CP-PACS Collaboration), Phys. Rev. D 64 (2001) 074510 [hep-lat/010302 8]. 

[4] Y. Iwasaki, K. Kanaya, S. KayaandT. Yoshie, Phys. Rev. Lett. 78(1997) 179 [hep-lat/960 902 2]. 

[5] O. Kaczmarek, F. Karsch, E. Laermann, C. Miao, S. Mukherjee, P. Petreczky, C. Schmidt, 
W. Soeldner, and W. Unger, arXiv: 1011 . 3130 [hep-lat] 

[6] M. Bochicchio, L. Maiani, G. Martinelli, G. Rossi, and M. Testa, Nucl. Phys. B 262 (1985) 331. 

[7] D. Toussaint, Phys. Rev. D 55 (1 997) 362 [hep - 1 at / 9 6 7 8 4 ] . 

[8] Y. Maezawa, N. Ukita, S. Aoki, S. Ejiri, T. Hatsuda, N. Ishii, and K. Kanaya (WHOT-QCD 
Collaboration), Phys. Rev. D 75, 074501 (2007) [hep-lat / 7 02 4]. 

[9] S. Ejiri, Y. Maezawa, N. Ukita, S. Aoki, T. Hatsuda, N. Ishii, K. Kanaya, and T. Umeda (WHOT-QCD 
Collaboration), Phys. Rev. D 82 (2010) 014508 [arXiv : 9 9 . 2 12 1]. 



7 



